Capillary Rise of Water in Hydrophilic Nanopores 
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We report on the capillary rise of water in three-dimensional networks of hydrophilic silica pores 
with 3.5 nm and 5nm mean radii respectively (porous Vycor monoliths). We find classical square 
root of time Lucas- Washburn laws for the imbibition dynamics over the entire capillary rise times 
of up to 16 h investigated. Provided we assume two preadsorbed, strongly bound layers of water 
molecules resting at the silica walls, which corresponds to a negative velocity slip length of -0.5 nm 
for water flow in silica nanopores, we can describe the filling process by a retained fluidity and 
capillarity of water in the pore center. This anticipated partitioning in two dynamic components 
reflects the structural-thermodynamic partitioning in strongly silica bound water layers and capillary 
condensed water in the pore center which is documented by sorption isotherm measurements. 
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Water transport in nanoscale environments plays a cru- 
cial role for phenomena ranging from clay swelling, frost 
heave, and oil recovery to colloidal stability, protein fold- 
ing, and transport in cells and tissues [Tl. The mobility 
of water has, therefore, been extensively studied by ex- 
periment and theory in such restricted geometries over 
the last three decades. 

These studies have revealed a remarkable fluidity of wa- 
ter down to nanometer and even subnanometer spatial 
conflnement \^ and have also demonstrated the validity 
of macroscopic capillarity concepts at the nanoscale and 
mesoscale However, how capillary forces along with 
the retained fluidity of water and other liquids contribute 
to the huge variety of phenomena, where self-propelled 
fluid transport is encountered in nanoscale geometries 
has been almost solely investigated so far by theory or 
for simple geometries as thin fllms [1]. 
Here, we report an experimental study on the sponta- 
neous imbibition, that is, the capillary rise (CR) of water 
in networks of hydrophilic silica pores with characteris- 
tic radii of 0(10 nm). Our study, which extends a for- 
mer more qualitative investigation of this phenomenology 
[5] , is aimed at a quantitative description of this process 
in terms of the driving capillary forces, the fluidity of 
water, and the complex flow topology encountered. We 
shall also demonstrate that an understanding of the CR 
dynamics requires a detailed knowledge of the thermo- 
dynamic state of water in this nanoscale environment, 
which we achieved by sorption isotherm measurements. 

As porous host we have chosen "thirsty" Vycor glass 
(Corning glass, code 7930) [IT, a virtually pure fused sil- 
ica glass permeated by a three-dimensional network of 
interconnected tortuous pores (see Fig. [1]). The exper- 
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FIG. 1: (Color online), (a) Schematic side view on the capil- 
lary rise of water in comparison with a raytracing illustration 
of spontaneous water imbibition in Vycor (b). A liquid col- 
umn has advanced up to the height h{t). A parabolic fluid 
velocity profile along with preadsorbed water layers beyond 
h(t) and a shaded resting boundary layer are sketched for the 
nanocapillary in panel (a). 



iments were performed with two types of Vycor signifl- 
cantly differing only in the mean pore radius (3.4 nm and 
4.9 nm respectively). We cut blocks of ~ 30mm height 
from the rods. Prior to using, we subjected them to 
a cleaning procedure with hydrogen peroxide and nitric 
acid followed by rinsing in deionized Millipore water and 
drying at 200 °C in vacuum for 2 days. Infrared spec- 
troscopy then reveals that the silica surfaces are covered 
by silanol groups with ^ 4 Si-OH/nm^ surface density 
[7], rendering them hydrophilic with a contact angle for 
water droplets close to [8J. At last, all facets but the 
bottom one were sealed. Water and nitrogen sorption ex- 
periments were performed with crunched Vycor samples 
at 4°C and -196 °C, respectively, using an all-metal gas 
handling setup. 

The mass uptake of the porous host as a function of im- 
bibition time m{t) was recorded gravimetrically with an 
electronic balance after immersing the monolith into a 
water bath thermostated to T = 24 °C. Right upon con- 
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tact of the bottom facet of the sample with the surface 
of the bulk liquid, we observed a jump in m(t) of about 
0.3 g. It can be attributed to the formation of a macro- 
scopic liquid meniscus at the perimeter of the matrix. 
Since we are interested in the water immersion into the 
inner pore space, we will present data sets after subtract- 
ing this mass jump. In principle, one has also to worry 
about buoyancy forces acting on the porous host and in- 
ertia effects in the initial states of the imbibition process 
[9] . However, both contributions are negligible in our ex- 
periment; a statement supported by the data presented 
below. 

As can be seen in Fig. |2^, we find a monotonically in- 
creasing m{t) up to i ^ 7h and 16 h, respectively. At that 
time, the advancing water front reaches the top of the 
monolith, as can be observed by the eyes due to a small 
contrast in transparency between the empty and filled 
parts of the sample. The time behavior of the imbibition 
of water into the porous monolith is well known and fol- 
lows a Lucas- Washburn law [TD] , which results from a 
competition of time invariant capillary forces driving the 
flow and an increasing viscous drag due to the increasing 
length of the paths connecting the advancing water front 
with the bulk reservoir. For macroscopic capillaries, such 
a simple scaling is only observable in the initial states of 
an imbibition process, eventually gravity comes into play 
leading to an exponential relaxation toward a final CR 
height characterized then by a balance of capillarity and 
gravity. By contrast, for liquids, wetting the inner sur- 
faces of capillaries with radii tl of a couple of nanometers, 
the driving capillary pressure, given by Laplace's law. 
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can be on the order of several hundred bars. In particu- 
lar, for water with its surface tension a of 72mN/m, Eq. 
[T] yields a pL ~ 300 bar for a pore with radius tl = 5nm. 
This corresponds to a hypothetical CR height of 3 km 
|llj . Thus, gravity forces should be negligible in our ex- 
periment and, accordingly, we find in accord with most 
recent microscopic simulations for liquid imbibition into 
nanopores ^ an excellent agreement with the Lucas- 
Washburn prediction over the entire t-range investigated 
(see the minute residuals shown in Fig. [2|d). 

More detailed quantitative insights regarding our CR 
experiment can be achieved by referring to the physics of 
liquid flow in porous media. According to Darcy's law, 
the volume flow rate V through a surface area ^ of a 
porous host which is subjected to a hydrostatic pressure 
Ap applied along the height h is given by 
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where rj is the shear viscosity of the liquid, and K refers 
to the hydraulic permeability of the porous matrix. 
For many porous systems with reasonable homogeneous 
cylindrical pore structures, among them also for Vycor, 
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FIG. 2: (Color online), (a) Normalized mass uptake m{t) of 
Vycor with two different pore radii (ro — 4.9 nm; triangles, 
ro — 3.4 nm: circles) due to water imbibition as a function 
of time t in comparison with Lucas- Washburn y^-fits (solid 
lines). Inset: imbibition speed Vi - tortuosity r map. The 
lines represent CR speeds determined with Eq. [5] from the 
measured m{t) rate as a function of the tortuosity r and for 
two different hydraulic pore radii (rh = ro and rh = ro — 2d, 
respectively, where d = 0.25 nm denotes the diameter of a 
water molecule) for both types of Vycor as indicated. The 
predicted Ui value for bulk water and its error margins are 
represented by horizontal lines, (b) Residuals Am between 
the measured m(t)-curve and the \/f-fits. 



it has been demonstrated that K can be well related to 
phcnomenological quantities typical of the matrix's mor- 
phology [Ellis], 



K = 
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where 0o is the volume porosity of the host, and r is its 
tortuosity characterizing the connectivity and meander- 
ing of the pores. Note that the hydrodynamic radius rh, 
which is the radius over which a parabolic flow proflle 
is established (see Fig. [TJi), does not necessarily have 
to agree with the mean pore radius tq because of either 
strongly adsorbed, immobile boundary layers, or due to 
velocity slippage at the pore walls [TH . In a CR ge- 
ometry, the liquid encounters the pressure drop along a 
time-dependent length h(t) (see Fig. [T|d), and Ap is solely 
determined by the Laplace pressure Ap — pi^ that acts at 
the advancing water menisci. Moreover, with the initial 
porosity t/fj of the sample upon CR start the permeated 
volume can be determined to be V{t) = (t)iAh{t). Ac- 
cordingly, Darcy's law transforms to a simple differential 
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equation, 
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Equation ^ is solved by a -y/t law for V, which by mul- 
tiplication with the mass density p of the liquid yields 
the aforementioned Lucas- Washburn law for the mass up- 
take. 
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As can be seen in Eq. [5) the ratio of surface tension to 
viscosity of pore-confined water, to which we refer in the 
following as imbibition speed Vi, is the property of the 
liquid which determines the dynamics of the CR. Thus, 
we can determine this ratio for confined water from our 
CR experiment provided we can assess the geometry 
factor G of our porous host. 

The parameters (j)o ^^nd rg are accessible by means of 
a gas sorption isotherms, which motivated us to record 
nitrogen and water vapor isotherms at T = —196 °C and 
T — 4°C, respectively. The latter is presented in Fig. 
|3] We plot the filling fraction /, that is, the number 
of water molecules adsorbed by the matrix normalized 
to the water amount necessary for its complete filling, 
versus the reduced vapor pressure p = Pc/po- The 
pressure po refers to the bulk vapor pressure of water at 
T = 4°C and Pc to the equilibration pressure after each 
adsorption or desorption step, respectively. The initial 
steep increase of / with p up to point A indicated in 
Fig. |3] is due to the formation of a first layer of water 
molecules which is strongly bound to the silica surface 
and hence exhibits a small p. Along the path A^B, a 
second layer grows. After the formation of the third layer 
in point C (/ ~ 0.3), we enter the hysteretic capillary 
condensation/evaporation regime [12], [(C)-(E)]. The 
entire pore space fills or empties via the formation or 
retreat of capillary bridges (see the inset in Fig. |3]). At 
point F, the matrix is completely filled and we measure 
Po typical of bulk droplets sitting outside the pores. The 
overall number of adsorbed molecules yields a porosity 
of (31 ±2)%. The mean pore diameters were determined 
from an analysis of the nitrogen sorption isotherms to 
(6.8 ± 0.2) nm and (9.8 ± 0.2) nm based on a mean-field 
model for capillary condensation proposed by Saam and 
Cole [16]. It is important to note that, according to the 
isotherms, the final relative humidity in our laboratory 
of ~ 30% {p = 0.3) will lead to at least one layer of 
strongly adsorbed water. This also means that the initial 
porosity 0i in our CR experiments is reduced by ~ 10%. 
Additionally, we have to assume tl to be reduced in 
comparison to vq by the thickness of the preadsorbed 
water film. 

The tortuosity could be inferred from experiments 
on the self-diffusion of liquids in porous Vycor [15) to 
r = 3.6 ± 0.4, in accordance with simulations of its 
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FIG. 3: (Color online). Water isotherms recorded at T = 
4°C for Vycor with 9.8 nm (triangles) and 6.8 nm (circles) 
mean pore diameter. Plotted is the filling fraction / versus the 
reduced vapor pressure p. Five characteristic points discussed 
in the text are indicated by solid symbols. 



pore morphology '17j. A value of r in the proximity 
of three seems reasonable if one considers that in an 
isotropic spongelike medium such as Vycor, the porosity 
can, in first approximation, be accounted for by three 
sets of parallel capillaries in the three spatial directions; 
whereas only 1/3 of these capillaries sustains the flow 
along the pressure drop Ap. A value larger than three 
reflects then the extended length of a meandering 
capillary beyond that of a straight one. 
We can now determine the CR speed Vi from our mass 
uptake measurements. In order to allow the reader a 
judgement on the influence of the uncertainties in r and 
Th on Vi, we plot the derived imbibition speeds as a 
function of t for selected assumed hydraulic radii rh (see 
the inset in Fig. |2|l). In principle, all of the resulting 
lines run in the vicinity of Vi = (79 ± 2) ™ calculated 
from bulk viscosity and surface tension of water at 
T = 24°C and indicated by the horizontal lines. More 
interestingly, for t = 3.6 ± 0.4, both experiments repro- 
duce the predicted bulk value for Vi {vi — (73 ± 40) ™ 
and Vi — (80 ± 33) ™ for Vycor with tq — 3.4 nm and 
ro = 4.9 nm, respectively) provided we assume two 
monolayers of water adjacent to the pore walls to be 
immobile (rh = rg — 2d, where d = 0.25 nm refers to 
the diameter of a water molecule). Note, for r\^ — ro, we 
systematically underestimate Vi in both cases. Thus, our 
experiments confirm former findings on the conserved 
fiuidity [5] and capillarity |3] of nanoconfined water, if 
we assume a velocity slip length; that is, the distance 
with respect to the fluid/solid interface, where the 
flow velocity vanishes, of — 0.5nm. This assumption is 
corroborated by recent molecular-dynamics studies on 
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the glassy structure of water boundary layers in Vycor 
and the expected existence of sticky boundary layers 
in Hagen-Poiseuille nanochannel flows for strong fluid- 
wall interactions [TS]. Also tip-surface measurements 
document a sudden increase of the viscosity by orders 
of magnitude in 0.5 nm proximity to hydrophilic glass 
surfaces jl9j. It also extends former experimental find- 
ings with respect to the validity of the no-slip boundary 
condition for water/silica interfaces, where this condition 
has been proven down to at least 10 nm from the surface 
[2U] . Moreover, the retained bulk behavior agrees with 
our observation of a correct relative scaling of measured 
and expected Vi values for water and a set of n-alkanes 
upon invasion in nanopores|5]. 

Our experiments yield additional rheological details. As 
mentioned above, Eq. [3] intrinsically assumes a parabolic 
flow velocity proflle across the pore cross-section, which 
implies a linear variation in the viscous shear rate, start- 
ing with in the pore center to a i-dependent maximum 
at Th, 7in oc For example, for the rh — 2.9 nm 

experiment, we estimate 7in to decrease from 7 x 10^ ^ 
after 1 s to 3 x 10^ - at the end of the CR. As we found no 

s 

i-dependent and, therefore, no 7-dependent deviations of 
m{t) from a single -v/t-fit, our measurements also testify 
the absence of any non-Newtonian behavior of water 
as well as an unchanged no-slip boundary condition 
imposed at rh. Despite the relatively large 7's probed 
here, the latter is not too surprising. The viscous forces 
of 0{rjCp^) can only overcome the strong water/silica 
interactions of 0{A/d) (Hamaker constant A ^ 10"^^ J) 
for 7 > 10^^ ^ [13 m] - significantly beyond the 7's 



probed here. 

It is worthwhile to note that water encounters a negative 
pressure upon CR which linearly decreases from — pl at 
the advancing menisci to atmospheric pressure at the 
bulk reservoir. Water's hydrogen bridge bond network 
is expected to be responsible for an increase in 77 and 
decrease in p under such large tensile pressures. Based 
on thermodynamic models for stretched water |22j . we 
estimated an ~ 3% Vi decrease due to this effect, which 
is, unfortunately, well below our error margins in u;. 
In summary, we have demonstrated classical Lucas- 
Washburn CR dynamics for water in networks of 
hydrophilic silica nanopores by rather simple gravi- 
metric imbibition experiments. They can be described 
in reasonable way by macroscopic hydrodynamics 
provided a sticky preadsorbed boundary layer of about 
two monolayers of water molecules is considered. We 
hope that the CR study presented here will stimulate 
further experiments on the interplay of fluidity and 
capillarity on the nanoscale not only for water but also 
for other complex liquids [23] . In particular, we envision 
experiments on carbon nanotubes with graphitic less 
water attractive walls, where single file fiow and velocity 
slippage may affect the CR dynamics markedly 
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